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The evolution of primordial black holes in a flat Friedmann universe with a massless scalar field is 
investigated in fully general relativistic numerical relativity. A primordial black hole is expected to 
form with a scale comparable to the cosmological apparent horizon, in which case it may go through 
an initial phase with significant accretion. However, if it is very close to the cosmological apparent 
horizon size, the accretion is suppressed due to general relativistic effects. In any case, it soon gets 
smaller than the cosmological horizon and thereafter it can be approximated as an isolated vacuum 
. . . solution with decaying mass accretion. In this situation the dynamical and inhomogeneous scalar 

'/^ , field is typically equivalent to a perfect fluid with a stiff equation of state p = p. The black hole 

■ mass never increases by more than a factor of two, despite recent claims that primordial black holes 

' might grow substantially through accreting quintessence. It is found that the gravitational memory 

(N ■ 

scenario, proposed for primordial black holes in Brans-Dicke and scalar-tensor theories of gravity, is 
highly unphysical. 

■ 

^ ' PACS numbers; 04.70.Bw, 97.60.Lf, 04.25.Dm, 95.35.+d 

T— I ; I. INTRODUCTION 

^ , The observed anisotropy and polarisation of the cosmic microwave background radiation give very precise informa- 
' tion about the cosmic history since last scattering. Using a simple physical model, we can determine cosmological 
' parameters with very high accuracy and also learn about structure formation and reionisation. Nevertheless, it is 
T-H ■ still difficult to obtain information about the universe before big bang nucleosynthesis without making many extra 
' assumptions. In this context, primordial black holes (PBHs) could be one of the most important fossils of the very 
early universe. Such black holes may have formed directly from primordial density perturbations and may con- 
tribute to the cosmological 7-ray background radiation, the cosmic ray flux and the dark matter density. This leads 
to important observational constraints on the number of PBHs 2] and hence on models of the very early universe. 

To understand these constraints, it is important to know whether PBHs can accrete enough to become much more 
massive than they were at formation. This topic has a long history. It was originally claimed, on the basis of a simple 
Q ■ Newtonian argument, that a PBH much smaller than the cosmological particle horizon at formation would not accrete 
' much at all, whereas one comparable to the horizon size would continue to grow at the same rate as the horizon until 
CZ3 i the end of the radiation-dominated era . Since one would expect a PBH to be of order the particle horizon size at 
formation in order to collapse against the pressure, this suggests that all PBHs might grow to the horizon mass at 
^ , the end of the radiation era, which is about 10^^ solar masses. Since there is no evidence for such enormous black 

holes in the universe today, the conclusion seemed to be that PBHs never formed. 
^ . , However, this conclusion was very suspect. If PBHs really could grow as fast as the particle horizon, there would 
' have to exist spherically symmetric self-similar solutions to Einstein's equations containing black holes in an exact 
Friedmann background. A study of solutions of this kind containing radiation (i.e. fluid with equation of state 
p — p/3) showed that there are no self-similar PBH solutions if the black hole is formed by purely local processes 
(i.e. if the background is exactly Friedmann beyond some radius). Self-similar solutions are possible only if the initial 
perturbation of the Friedmann background extends to infinity |j] . Since the PBH must soon become much smaller 
than the cosmological horizon, at which point the Newtonian argument should be applicable, this suggests that it will 
only increase its mass by a small amount. 

This conclusion was subsequently extended to more general fiuids, with an equation of state of the formp = kp \ 
Although there was an initial claim that self-similar growth might be possible in the special case of a stiff fiuid 
{p — p) 1 this claim was subsequently challenged 8| . It was found to be true only in rather contrived circumstances 
in which the stiff fluid is converted into radiation at the black hole's event horizon. Therefore the conclusion that there 
are no self-similar solutions containing PBHs formed by purely local processes seems to be true very generally. See 
for a review of self-similar solutions. This conclusion is also confirmed by numerical calculations. If the formation and 
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evolution of a PBH in a general fluid universe with a local perturbation is simulated numerically, without assuming 
self-similarity, it is found that the PBH soon becomes much smaller than the cosmological horizon and this excludes 

PBH formation from density perturbations has been investigated in the 
but again with no evidence for self-similar growth. 
Hitherto studies of PBHs have mainly focused on perfect fluid universes with equation of state p = kp. However, 
it is also natural to consider a universe whose density is dominated by a scalar field. For example, in the chaotic 
inflation scenario it is postulated that there is a pre-inflationary stage in which the scalar field moves randomly in 
space and time. In the pre-heating inflationary scenario, it is also natural to consider a scalar-field-dominated era Jl'^. 
More recently, the study of PBHs in the quintessence scenario has attracted attention . It is therefore important 
to examine whether the conclusion that PBH accretion is small also applies in this case. 

If the scalar field is massless and there is no scalar potential, then it is well known that it is equivalent to a stiff 
fluid providing the gradient of the scalar field remains timelike, as usually applies. The fact that there is no self-similar 
solution in the stiff case therefore suggests that accretion is also limited in the scalar field case. Despite this, the 
original Newtonian argument for self-similar growth has recently been applied in the quintessence scenario to argue 
that PBHs could grow enough to form the supermassive black holes found in galactic nuclei A generalization 

of this analysis _21j . not necessarily involving self-similarity but still based on the Newtonian analysis, has also 
claimed there could be appreciable quintessence accretion in some circumstances. Both these analyses contravene the 
conclusions of the earlier work discussed above. However, since the earlier work did not strictly include the scalar 
field case and did not allow for a scalar potential, a more careful analysis is required before concluding that these 
analyses are erroneous. 

Scalar fields are also relevant in Brans-Dicke and scalar-tensor theories of gravity, where the gravitational "constant" 
varies in space and time. This is because, if there is a sin gle g ravitational scalar field, such theories can be transformed 
into the usual Einstein gravity with a single scalar field |22l |. Brans-Dicke and scalar-tensor theories are particularly 
relevant to PBHs since the black holes may form when G was very different from today. Indeed the "gravitational 
memory" scenario has been proposed in which the value of G within the black hole is assumed to be preserved 
as the cosmological background value evolves. The observational constraints on PBHs depend strongly on whether 
or not one has gravitational memory 24] . It is not clear whether this applies but, if it does, it should be due to the 
properties of the black hole event horizon rather than those of the matter fields involved. 

It is interesting that the problem of gravitational memory is also closely related to the problem of accretion |25j | . 
This is because it turns out that appreciable accretion of the scalar field energy is required for gravitational memory 
to be preserved and so the issue of whether there is a self-similar solution again becomes relevant. Indeed Carr and 
Goymer argued that there is unlikely to be gravitational memory on the grounds that there is no self-similar solution 
in the stiff fluid case |25j . Evidence against t he g ravitational memory scenario, at least in particular situations, has 
also been obtained using other arguments |26l l27j . 

As a first step to investigating these questions numerically, in this paper we consider a universe containing a massless 
scalar field but no matter. To implement the simulations, we use a double-null formulation of the Einstein equations 
based on the work by Hamade and Stewart [2^. This has been shown to be a very powerful tool for investigating 
critical collapse |28j[ and the internal structure of black holes |29| . This study represents an improvement on earlier 
work in which the gravitational effect of the scalar field was neglected, since it uses the full field equations. It 
should also be seen in conjunction with two other papers [30l Isij : we consider general constraints on the size of a 
PBH at formation in the first and some effects associated with very large PBHs in the second. 

The plan of the paper is as follows. In Section II, we describe the double-null formulation and present the basic 
equations for the scalar field case. In Section HI, we discuss how we set up the initial data for the simulations. Section 
IV presents our results, with particular emphasis on the evolution of the black hole event horizon and the spatial 
profile of the scalar field. In Section V, we discuss the implications of these results for the accretion and gravitational 
memory issues, also considering the applicability of the Newtonian accretion formula. In Section VI, we draw some 
general conclusions. Details of the numerical code and particular exact solutions are described in Appendices. We 
adopt units in which G = c = 1 and the abstract index notation of reference !32| . 

II. DOUBLE-NULL FORMULATION OF THE EINSTEIN EQUATIONS 

We consider a massless scalar field in general relativity, for which the stress-energy tensor is 

Tafc - *,a*,b - ^5a6*''*,c. (2.1) 
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The Einstein equations are 




(2.2) 



and the equation of motion for the scalar field is 



□ * = * 



;a 



;a 



= 0. 



(2.3) 



We also focus on a spherically symmetric system, for which the line element can be written in the form 



ds^ = -a^(u, v)dudv + r^(M, v){dO'^ + sin^ 6d(j)^), 



(2.4) 



where u and v are advanced and retarded time coordinates, respectively, and r is the "area radius" (the proper area 
of the sphere of constant r being Airr^). Eqs. (|2.2I) and (|2.3I) then imply that we have 14 first-order partial differential 
equations and two auxiliary equations (see Section 2 of .28] ) . By adopting a double- null coordinate choice, we can 
simulate regions outside the cosmological apparent horizon and inside the black hole apparent horizon simultaneously. 

In spherically symmetric spacetimes the existence and position of apparent horizons can be inferred from the form 
of the Hawking mass. This is a well-behaved quasi-local mass, which can be written as 



A region is trapped if < and r^^, < 0, while it is antitrapped if r^u > and > 0. The black hole and 
cosmological apparent horizons are defined as marginally trapped and anti-trapped surfaces, respectively. Providing 
the black hole horizon is within the cosmological horizon, these correspond to the conditions r^y = and r „ = 0, 
respectively. Thus the relation r — 2m is satisfied on both apparent horizons because r.,,r.7 i = there. We confine 
attention to this situation in the present paper. However, as discussed in a separate paper |3l|, in some circumstances 
the black hole horizon can be outside the cosmological horizon and the situation is then more complicated. In this 
case, we can still use two different marginal surfaces on which r^^ = and = 0, respectively However, these 
are no longer everywhere identified with the black hole and cosmological apparent horizons. 

Note that the black hole and cosmological apparent horizons are distinct from the black hole event horizon and 
cosmological particle horizon, which are always null. For numerical purposes, the apparent horizons are easier to find 
than the event and particle horizons because the complete spacetime, including the initial singularity, is needed to 
identify the latter, whereas the former can be found from the quasi-local properties of spacetime alone. However, we 
will also discuss the behaviour of the black hole event horizon in Section Hvl 

We first consider the flat Friedmann model, for which a^(u,w), r{u,v) and m{u,v) are given explicitly in terms of 
the standard double null coordinates in Appendix^ The condition r = 2m implies that the cosmological apparent 
horizon is 3u -|- = and this implies that its radius is just the Hubble scale H^^ in a flat Friedmann universe, as 
described in '3^. On the other hand, the cosmological particle horizon is clearly u = 0, corresponding to a photon 
propagating outwards from the big bang. This implies that the apparent horizon is always outside the particle horizon 
and the conformal diagram of the spacetime is as indicated in Fig. ^ This also shows the initial (big bang) spacelike 
singularity. If one perturbs such a spacetime but without introducing a black hole, the conformal diagram remains the 
same but the trajectories of constant space and time coordinate change. If one has a black hole embedded in an exact 
or asymptotically flat Friedmann model, the conformal diagram will change to the form indicated in Fig. |21 There 
are now also a black hole event horizon and apparent horizon, although one cannot give explicit analytic expressions 
for these, and a final (black hole) spacelike singularity. 

If we define time and space coordinates, T and X , by 




(2.5) 



Using the above equation and Eq. (2.6) of [23, we can derive the following useful relations: 

8^r2r,„(*.„)2 
m.u = 5 — '■ , 

8^r2r„(*^,)2 



(2.6) 



(2.7) 



u 



V 



T-X, 
T + X, 



(2.8) 
(2.9) 



we can put the 2-dimensional part of the metric tensor into conformally flat form. As in the usual 3-1-1 approach, we 
define the energy density pu and momentum density J measured by an observer moving normal to the T = const 
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spacelike hypersurface by 

PH = T„fcf'^f^ (2.10) 

J = -Tabf''X\ (2.11) 

where and X'^ are unit vectors parallel to T° and X^, respectively. 

It should be stressed that the above description is observer-dependent and one can adopt another point of view. 
Providing the scalar field is vorticity-free and has a timelike gradient, it is equivalent to a perfect fluid with a "stiff" 
equation of state p — p 19]. The stress-energy tensor is then 

T"'' = psig'''' + 2U''U''), (2.12) 

where the energy density ps and the 4-velocity C/" of the stiff fluid are 

PS = -^*,a*^", (2.13) 
C/a = . (2.14) 

It should be noted that the 4-velocity of the stiff fluid is parallel to the (timelike) gradient of the scalar field. We can 
define an observer-independent velocity function by 

F=^=Vl7", (2.15) 

this being the rate of increase of r per unit proper time t along the worldlinc of the equivalent stiff fluid element. 
Appendix gives expressions for these physical quantities in terms of the quantities calculated numerically. 



III. INITIAL DATA FOR PBHS 



A. Structure of initial data 



The initial data are prescribed on the outgoing null surface u — uq and the ingoing null surface v ~ vq. The region 
of calculation is the diamond [uo,ui] x [fo,^^i] shown in Fig.O which is also equivalent to the diamond in Fig.|2] We 
have three independent functions on the two null surfaces: a^, ^ and r. Two of them can be chosen freely and the 
other one is determined by the initial value equations on the null surfaces. It is convenient to choose 

a^(uo, w), a^(u, wq), *("o, «), *(«, "o) (3.1) 

as the free initial data and to regard 

r{uo,v),r{u,vo) (3.2) 

as being determined by the initial value equations. We can regard ^'(uo,w) and ^{u,vo) as the physical degrees of 
freedom in the initial data, while the choice for a^{uo, v) and a^(it, vq) fixes the gauge. 

It is usually assumed that the perturbation from which the PBH forms is local, so that the universe is exactly 
Friedmann at sufficiently large distances. Since information propagates at the speed of light, the matching between 
the perturbed and flat Friedmann regions always corresponds to some outgoing null ray u = Um = const and such a 
ray is shown in Fig. O This means that the region outside the matching surface u — is always described by the 
flat Friedmann solution. On the outgoing null surface u = uq, we therefore assume that the initial data are given by 
the flat Friedmann solution for vq < v < vi. On the ingoing null surface v = vq, we assume they are given by the flat 
Friedmann solution for uq < u < w,„ but by some perturbation of it for < u < ui, where an appropriate matching 
is implemented at w = Um- As described in the next section, the perturbed region contains a Schwarzschild black 
hole. However, since this is a vacuum solution, it only applies instantaneously at v — vq because the inflowing matter 
will fill the vacuum immediately. 

If the perturbation is assumed to be local, we can derive various upper limits on the black hole size |30l | . These are 
generally within the cosmological particle horizon, which would then correspond to the line u = in Fig.|21 However, 
a perturbation can be larger than the cosmological particle horizon in some circumstances - indeed most of the cases 
we consider in this paper have this feature. The particle horizon is then within the perturbed region and no longer 
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given by It = 0. In some circumstances, e.g. for quantum fluctuations resulting from inflation, the perturbation may 
extend indefinitely, or at least well beyond the usual particle horizon. Generally the only upper limit on the size of 
the perturbation comes from the requirement that the the perturbed region must be part of our universe rather a 
separate closed universe. The condition for this has been derived precisely for the situation in which the collapsing 
region is homogeneous and the equation of state is p = kp this case, the closed-universe scale is always well 

beyond the cosmological apparent horizon size. Even if the perturbed region is within our universe, the solution may 
exhibit anomalous features if it is too large. These anomalies are described elsewhere ^31] but we do not choose initial 
data which lead them to arise in this paper. 



B. Setting up the initial data 



In the flat Friedmann region, we choose the coordinate system given in Appendix lB II We also impose flat Friedmann 
initial data for and ^ on the initial outgoing null surface u = uq with vq < v < vi. This is also given in Appendix lB II 
On the initial ingoing null surface v = vq, we choose flat Friedmann data for for uq < u < ui. For the function 
we use the same data on the initial ingoing null surface for uq < u < Um, but 'J = const for Um + Am < u < ui. 
This is equivalent to Schwarzschild data in coordinates penetrating the black hole (cf. Appendix IB 2(1 . The sudden 
transition from flat Friedmann data to Schwarzschild data results in a discontinuity at u = u^, which reduces the 
Hence we smooth the transition with some smoothing length Am; we use a quadratic function 
f Am, so that and '^^u are continuous. More precisely, we impose the following initial data for 



numerical accuracy, 
between Mm and Mm 
and s = Vivr*: 



2 ( u + Vo 
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(3.3) 

(u < Mm) 

So (Wm < W < Mm + Am) , (3.4) 
(m > Mm + Am) 

(3.5) 
(3.6) 



on the initial outgoing null surface m = Mq. Here C and sq are constants and, without loss of generality, we can choose 
C = 1 and So = 0. Note that one has a Schwarzschild vacuum for u > Mm + Am. However, as noted above, this 
situation only applies instantaneously at m = mo, since the vacuum will be filled with the scalar field at any later or 
earlier time. Nevertheless, one expects the event horizon to evolve smoothly, so it will still be described by a single 
null line, as illustrated qualitatively in Fig. El 

To summarise, two parameters describe the initial conditions: the location of the matching surface and the smooth- 
ing length. The region outside the matching surface is the exact flat Friedmann solution, so the first parameter 
determines the mass of the perturbed region. The second parameter determines the mass of the black hole compared 
to the perturbed mass. A small smoothing length means a narrow boundary between the exterior Friedmann region 
and the interior vacuum region. In this case, the mass inside the vacuum region is nearly the same as that within the 
matching surface and so a large mass is concentrated at the central singularity. A large smoothing length means a 
wide boundary, so the mass inside the vacuum region is much smaller than this and so a small mass is concentrated at 
the central singularity. For vanishing smoothing length, the black hole mass would be the same as that at the match- 
ing surface, so Um would also be the black hole event horizon, although it is difhcult to implement this numerically 
because of the strong discontinuity. 



IV. RESULTS 

The numerical code is based on Hamade and Stewart although it is slightly modified so as to be accurate to 
second-order (see Appendix . As indicated in Fig. |31 the calculated region is the diamond contained by mo, mo, ui 
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and vi. The initial data are prescribed on the ingoing null surface vq = 1 and the outgoing null surface uq = —0.5. 
The values Ug and Vq fix the units. For example, the background value for the Hubble parameter at {uo,Vo) is 4 and 
the Hubble length becomes 0.25 from Eq. ljB8(l . The value of ui is always taken to be 1.1 but vi varies between 2 and 
4. On the initial ingoing null surface, we make the matching at three values of Um in the range —0.5 to 0. We also 
use small and large smoothing lengths, corresponding to Au — 0.02 and Au = 0.5, respectively. 

As shown in Appendix IB II the cosmological apparent and particle horizons are given by ucah = — f/3 and 
WCPH = 0, respectively, in the flat Friedmann model. As time proceeds, the outgoing null rays will become ever more 
sensitive to r near the black hole event horizon, so the calculation is stopped at the value of v at which the outgoing 
null rays near the black hole event horizon become too coarse to resolve. This then specifies the value Vi in Fig. O 
As shown in Appendix IB II the radius of the unperturbed cosmological apparent horizon and the mass within it are 
given by 

mcAH=rcAH/2 = Ci;^/V(3\/3). (4.1) 

However, this equation must be modified if the perturbed region extends beyond the usual cosmological apparent 
horizon and this applies in one of the cases considered. 

The model parameters and initial black hole to cosmological horizon mass ratios are summarised in Table ^ We 
present results for four models. For Model A, a sharp matching is made at the cosmological apparent horizon. For 
Model B, a sharp matching is made at the cosmological particle horizon, which is well inside the cosmological apparent 
horizon. For Model C, a smooth matching is made at the cosmological apparent horizon. For Model D, a smooth 
matching is made well outside the cosmological apparent horizon, although the black hole itself is inside it. In a 
separate paper |3lj| . we consider models in which a sharp matching is made well outside the cosmological apparent 
horizon. These models have very different qualitative features, even the conformal diagram being modified, which is 
why we consider them separately. 

A. Black hole horizons 

Figure 0] shows the locations of the black hole event horizon and apparent horizon and the cosmological apparent 
horizon in the {u, v) plane for the four models. Each figure can be identified with the diamond in Fig.|31 tilted through 
45 degrees. In all cases, it is seen that the black hole apparent horizon is initially well within the event horizon but 
approaches it as time proceeds. In Fig. 2Ib), the location of the (unperturbed) cosmological particle horizon of the 
exact flat Friedmann solution is also plotted. In other cases, the cosmological particle horizon is perturbed, so we do 
not plot it. The cosmological apparent horizon is unperturbed except in case D. 

The radii and masses of the black hole apparent horizon and event horizon for Models A-D are shown, together with 
those for the cosmological horizons, in Figs.[Sland|^ The initial and final values and their ratios are also summarised 
in Tables HTl and IIIII All three horizons approximately coincide at ?; = 1 for Model A, both in radius and mass, as seen 
in Figs. Ela) and Efa). Figs. Eland El show that the qualitative features of the evolution for Models A-D are very 
similar. The radius of the cosmological apparent horizon of the unperturbed flat Friedmann solution grows like u"^/^, 
while the black hole apparent horizon and event horizon converge and grow much more slowly. As a consequence, 
the black hole horizons soon get much smaller than the cosmological apparent horizon. This is also the case for the 
corresponding masses. 

Since the black hole event horizon is given by the curve u — mbheh — const, Eq. H2.7|l gives the mass accretion 
rate of the black hole. This also serves as a consistency check for numerical accuracy. The result for Models A-D is 
shown in Fig.[7| where both sides of Eq. (|2.7|l are plotted. We find that the curves are almost indistinguishable, which 
implies very good numerical accuracy. For Model A, the accretion rate starts very small, then increases, reaches a 
maximum of 0.0691 at w ~ 1.75 and then decreases to a very small value. Recall that the units are such that the 
background Hubble parameter is 4 and the background Hubble length is 0.25 at uq and vq- For Models B and C, the 
accretion rate starts with its maximum value and then monotonically decreases. The features for Model D are similar 
to those for Model A. 

An important qualitative feature is that accretion is suppressed for a PBH nearly as large as the cosmological 
apparent horizon. (In Model D the black hole event horizon happens to be close to the cosmological apparent 
horizon.) This is because Eq. (|2.7|l implies that the mass accretion rate is proportional to (— r_„), i.e. the ingoing 
null contraction, and this vanishes at the cosmological apparent horizon. This feature is a purely general relativistic 
effect and there is no analogue in Newtonian gravity. It can be understood physically as arising because the scalar 
field must be partaking in the cosmic expansion at sufficiently large distances, rather than accreting onto the black 
hole. Since the sound speed is the speed of light for a scalar field, the cosmological apparent horizon is likely to be 
the transition point. A more rigorous explanation of this effect is presented in Section V, where it is shown to apply 
for more general fluids. 
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We conclude that, even if a PBH starts off as large as the cosmological apparent horizon, it soon gets considerably 
smaller than it. Whatever the initial size of the PBH, accretion can be large only for about a Hubble time, after 
which it becomes ineffective. As seen in Tables IhI and IHll the mass increase is at most 79 % for our models. 

B. Geometry and scalar field 

Since the qualitative features of the numerical results are very similar for Models A-D, in this subsection we 
concentrate on Model A. Figure |S1 shows the evolution of the area radius r in terms of u and v. We can see from 
Fig. IHIb) that there is a threshold value of u: outgoing null rays with smaller u go to infinity, while those with larger 
u go to the singularity at r = 0. Therefore, this threshold value can be identified with the black hole event horizon 
WBHEH — 0.702. Since r^^ approaches zero along the black hole event horizon, all ingoing null rays seem to cross it at 
almost the same radius in Fig.|S{a). In fact, the integration is continued well beyond v = vi to locate the black hole 
event horizon accurately. Figure Hlb) shows that the radius of the black hole event horizon reaches 0.694 dX v — A 
and this might be regarded as its asymptotic value. 

Figure El shows the evolution of 2m /r in terms of u and v. There is an apparent horizon where this is unity. The 
signs of both ingoing and outgoing null expansions are the same when this is larger than unity. A region is untrapped 
when it is smaller than unity. As seen in Fig.|^a), 2m /r is almost exactly unity in the perturbed region at the initial 
ingoing null surface v = = \. As u increases, there then appears a black hole apparent horizon, beyond which the 
region is trapped. There is also a cosmological apparent horizon at u = —1/3 and v = vo = I, although it goes outside 
the calculated region during the evolution. As time proceeds, the value of 2m /r in the untrapped region between the 
black hole and cosmological apparent horizons becomes smaller and the curve near the black hole apparent horizon 
becomes steeper. As seen in Fig. |^b), there is no black hole apparent horizon before the event horizon forms at 
u = ubheh — 0.702. At u = mbheh, 2m/r approaches unity from below as v increases. Soon after u = urheh, a 
black hole apparent horizon appears. 

Figure [TUI shows the evolution of the scalar field for Model A. For clarity, it is plotted in terms of both {v,r) 
and (r, r) coordinates. In either case there is no coordinate singularity at the black hole horizons as the ingoing 
Eddington-Finkelstein coordinate is well-behaved for the Schwarzschild black hole. In the {v, r) diagram of Fig. lior a). 
the broken and solid light curves give the results of the simulations and the evolution of the flat Friedmann solution, 
respectively. Because of the initial perturbation, the value of the scalar field in the perturbed region starts off much 
smaller than that of the flat Friedmann solution. Thereafter, the perturbed scalar field tends to evolve as in the flat 
Friedmann solution. The position of the event horizon is shown by the solid heavy curve. Since the radius of the black 
hole event horizon is 0.694 at w = 4, the region inside the black hole event horizon is also calculated accurately. The 
evolution of the scalar field is well described by the flat Friedmann evolution with a small perturbation both outside 
and inside the black hole event horizon. The scalar field is smooth at the black hole event horizon. Fig. Iior b) shows 
the equivalent results in the (T, r) diagram. It should be stressed, however, that there is no unique choice of spatial 
hypersurface and the profile of ^' would be different with another choice. 

Figure ITTI shows the energy density pn and momentum density J measured by the observer moving normal to the 
T = const spacelike hypersurface for Model A. Since this normal observer coincides with the comoving observer in 
the flat Friedmann universe, J vanishes in the flat Friedmann region for u < (see Appendix IB In terms of the 
energy density pu, there is an underdense region inside the flat Friedmann region. Inside this underdense region, there 
is an overdense region around the black hole event horizon, where the energy density increases with time. In terms 
of the momentum density J, there is considerable energy infall around the black hole event horizon. | J| increases 
significantly near the black hole event horizon, although we should note that this is an observer-dependent view. 

Figure El shows the observer- independent quantities ps and V. These are the energy density and velocity of the 
equivalent stiff fluid. The spikes correspond to the matching coordinate Um- It is interesting that ps is everywhere 
positive. This is a non-trivial result since ps could in principle be negative if the spatial gradients were large enough. 
This is related to the fact, as illustrated in Fig. I12r a'l. that the gradient of the scalar field is always timelike in this 
simulation. (Otherwise it would not be equivalent to a stiff fluid.) Also, in terms of the density pg, we can see an 
underdense region just inside the flat Friedmann region. Within this underdense region, the energy density increases 
as one moves inwards. However, the energy density decreases nearly homogeneously for v > 2.5. In terms of the 
velocity V, the region near the black hole is infalling, as seen in Fig. I12r b). but the infall velocity profile near the 
black hole event horizon is almost stationary. We conclude that the growth of pn and | J| around the black hole event 
horizon is completely due to the Lorentz factor. Physically, the energy density ps is slowly decreasing and the infall 
velocity profile V is almost stationary near the black hole event horizon. These results suggest that the PBH can 
be regarded as almost isolated, with a relatively small mass accretion for w > 2. This is consistent with the mass 
accretion rate seen in Fig.[7{a). Although there are apparent discontinuities at large radii in Figs. Illland ll2l this is 
only due to the narrow boundary and all quantities are continuous everywhere. 
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V. DISCUSSION 



A. Self-similar growtli and tlie Newtonian accretion formula 

Various arguments 0, ^ HE HOl lead to an estimate for the mass accretion rate of a PBH in a cosmological 
background of the form 



dt 



AnapVsrj^ (5.1) 



where p is the background density, is the sound speed, ta is the accretion radius of the hole and a is a dimensionless 
constant which depends on the equation of state {p = kp) and other features of the accretion flow. None of these 
arguments is precisely correct, since they all involve approximations which make them in some respects Newtonian. In 
particular, the value of a depends upon whether one allows for such effects as the relativistic sound speed, relativistic 
pressure and relativistic focussing. Also a stationary flow is usually assumed and the background cosmological 
expansion is neglected. Nevertheless it might be hoped that these inaccuracies could be absorbed into the constant 
a. 

If one assumes that the appropriate density to use in Eq. (|5.1|) is the background cosmological density ^ l/(Gt^), 
the accretion rate becomes 

~dr - 'W' ^^-^^ 

where (3 is another constant with units C^/G. This has the following solution: 

fit 

rriBH = — T- ^ (5.3) 

where Mf is the mass of the black hole when it forms at time tf.li the hole is much smaller than the particle horizon 
at formation {Mf ^ Ptf), this implies that there is very little growth. However, if Mf is chosen so that the term in 
brackets in Eq. (|5.3(l is zero, it suggests the possibility of self-similar growth with ttibh oc t. We assume the constant 
/3 is 6, since that is what the analysis of reference |8(l | implies when fc = 1. If the Newtonian formula were correct, 
the ratio would be order of unity. See for a full discussion of the solutions. 

In discussing the accretion rate for a scalar field, which is relevant here, one possibility is to regard it as a stiff 
fluid. In this case, one can apply the above analysis with /3 having the value appropriate for k ~ I and this is the 
approach used by Carr and Goymer 2^. Bean and Magueijo 20] also assume the accretion is given by Eq. H5.1|l but 
identify p with the kinetic energy of the scalar field (1/2)^*^. However, since this is assumed to falls off as they 
obtain an equation of the same form as Carr and Goymer, though with a different constant since they allow for a 
scalar potential. Custodio and Horvath allow for a broader range of behaviours for 4', leading to a generalization 
of Ea.f l5.3|l . All of these analyses suggest that the black hole could grow self-similar ly with a fine-tuning of the initial 
mass. However, as discussed in the Introduction, this conclusion conflicts with the analytic demonstration that there 
is no self-similar solution which contains a black hole attached to an exact flat Friedmann exterior solution 0, 0, S ■ 

The prediction of a large accretion rate also conflicts with the results of our simulations. Since the cosmological 
time is given hy t = (2/3)[(u -I- ■i;)/2]'^/^ for the flat Friedmann solution (see Appendix IB ip . the self-similar growth of 
a black hole would imply that its mass grew like {u + v)^^^. The black hole mass accretion rate dmBaEH/dv should 
therefore be proportional to {u + vY^^. However, the results of our simulations show that the accretion rate decreases 
after reaching a maximum. We conclude that there is no evidence for a trend towards self-similar evolution, which 
contradicts the argument that a black hole can grow as fast as the Universe through accreting a quintessence 
field. It also provides a counterexample to the proposal that a spherically symmetric spacetime always evolves towards 
self-similarity 0. 

It order to determine the accuracy of the Newtonian accretion formula more precisely, we now derive an exact 
relativistic formula for the accretion rate. Because the event horizon is given by an outgoing null surface u = const, 
the time variation of the black hole mass is governed by Eq. (|2.7|l . Also Eqs. (|A2|) and ljA3|) imply 

where pn and J are the energy density and momentum density. Combining this with Eq. (|2.7|l gives 

m^y = -Snr^r^uipH + J)- (5.5) 
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If we take the time coordinate to be 



t = h{u + w), 



(5.6) 



where h is an arbitrary function of u + w, then we have 

dv 



dm 



m . 



dt 



dr 
It 



du 



Combining this with Ea.r i5.5|l gives an expression for {dm/ dt) u=const- In particular, since the black hole event horizon 
has u =const, the accretion rate can be expressed as 



dm 



BHEH 



dt 



dr 
'dt 



(5.7) 



This accretion equation is exact and might be compared with the Newtonian prediction given by Eq. (|5.1(l . We see 
that p is replaced by pn + while c is replaced by the value of dr / dt along the path with constant v. 

Eq. H5.5|l also explains why the mass accretion rate starts low for models A and D. Since the black hole event 
horizon is inside but very close to the cosmological apparent horizon in these cases, is negative but very close to 
zero. After a while, however, falls well below zero and then the mass accretion increases. The suppression of the 
accretion by the factor on the right-hand side of Eq. (|5.5|l can be understood in a more general context from the 
counterparts of Eqs. (|2.6I) and H2.7|l for a general spherically symmetric spacetime [s^ils^- Using the present notation, 
the equations become 



^(T r -T r ) 

(2^ 



(5.8) 
(5.9) 



The combination of the two terms in parentheses on the right-hand side of Eq. (|5.9|l . related to outgoing and ingoing 
null expansions respectively, determines the time variation of the black hole mass. In the case of a massless scalar 
field, the situation is simphfied because T-av — , T„„ = (*,«)'^ > and Ty^ = {'^,v)'^ > 0- As a result, we can 
immediately conclude that the black hole accretion vanishes when its event horizon coincides with the cosmological 
apparent horizon, since = there. 

It can be easily proved that the black hole mass is non-decreasing if > 0, < and if the dominant energy 
condition holds on the event horizon (see Proposition 5 of |33|V However, if these assumptions are not satisfied, 
the mass variation of the black hole event horizon is non-trivial. In particular, one finds that the black hole mass 
is decreasing when it is outside the cosmological apparent horizon, i.e., r,u > 0. This possibility is studied in detail 
elsewhere [sj . It should be noted that Hayward [SJ, 13 HE HEl introduced the future (past) outer trapping horizon 
as a more general and useful concept than the black hole (white hole) event horizon. When the black hole (white 
hole) horizon is defined in this way, he also proved the monotonicity of the area and quasi-local mass in the spherically 
symmetric situation. 



B. Gravitational memory? 

The Brans-Dicke theory with empty stress-energy tensor can be transformed into the Einstein theory with a massless 
scalar field by a conformal transformation even when the Brans-Dicke parameter lo is of order unity . This is also 
the case in scalar-tensor theories of gravity if there is a single massless gravitational scalar field. The transformed 
frame with the Einstein-Hilbert action is called the Einstein frame, in contrast to the original physical frame. The 
gravitational constant varies in space and time in the physical frame and is actually a function of the massless scalar 
field \E' in the Einstein frame. Since the conformal transformation does not affect the causal structure of the spacetime, 
the present simulation can be regarded as probing the evolution of PBHs in Brans-Dicke and scalar-tensor theories 
with empty stress-energy tensor. 

Our results elucidate whether or not gravitational memory is physically reasonable p5l | . We have seen that the scalar 
field evolution near and even inside the black hole event horizon follows the cosmological evolution. Moreover, the 
gradient of the scalar field is always timelike. Therefore, we conclude that the evolution of the physical gravitational 
constant near a PBH in Brans-Dicke and scalar-tensor theories essentially follows its asymptotic cosmological evolution. 
The fact that the gradient of the scalar field is timelike during the evolution also ensures the existence of a spacetime 
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foliation with spacelike hypersurfaces, in which the physical gravitational constant is spatially constant. Therefore, 
we conclude that the gravitational memory scenario is physically unrealistic even for PBHs whose size is comparable 
with the cosmological horizon. 

Using a perturbative test field analysis, Jacobson p6| showed that black holes do not exhibit gravitational memory 
in scalar-tensor cosmology if they are much smaller than the cosmological horizon. Numerical studies for a test scalar 
field in the PBH dust solution by Harada, Goymer and Carr |23] showed that there is no gravitational memory even 
for black holes whose size is comparable to the cosmological horizon scale, so long as the back reaction of the scalar 
field on the metric can be neglected. In the present work we have extended this result by considering the case where 
the black hole size is comparable to the cosmological horizon scale and the gravitational effect of the scalar field is 
included. Once the black hole horizon gets much smaller than the cosmological horizon scale, the approximation 
made by Jacobson will apply. Taken together, these results suggest that the gravitational memory scenario is a highly 
unphysical, at least in Brans-Dicke and scalar-tensor theories. 



VI. CONCLUSION 



We have investigated the evolution of a PBH in a fiat Friedmann universe with a massless scalar field by numerically 
integrating the Einstein field equations using the double-null formulation. We have considered models in which a 
Schwarzschild interior and a Friedmann exterior are initially matched at a finite radius with some smoothing length. 
For PBHs which are initially the same size as or smaller than the cosmological apparent horizon, the black hole event 
horizon soon becomes smaller than the cosmological horizon scale. The black hole apparent horizon forms inside the 
black hole event horizon but asymptotes towards it. 

The scalar field evolution qualitatively follows the cosmological evolution close to the black hole event horizon and 
even inside it. The field maintains a timelike gradient and is therefore always equivalent to a stiff fiuid. In terms of 
the stiff fluid description, the energy density in the perturbed region decreases homogeneously in time and the infall 
velocity profile near the black hole becomes almost stationary after a short while. 

Our results show that, soon after the PBH enters the cosmological apparent horizon, it becomes isolated from the 
cosmological expansion, with only a very small amount of mass accretion from the non- vacuum exterior. The accretion 
can be significant at first but soon decreases and becomes insignificant. In particular, our simulations exhibit no self- 
similar growth of the PBH. This is consistent with earlier work, which proved the non-existence of self-similar PBH 
solutions for a stiff fluid with an exact flat Friedmann exterior. The present result also indicates that the gravitational 
memory scenario for PBHs is unphysical in Brans-Dicke and scalar-tensor theories of gravity. 
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APPENDIX A: PHYSICAL QUANTITIES AND NUMERICAL CODE 



Following Hamade and Stewart we use variables, s. A, c, d, /, g, p, q as well as auxiliary variables A and /i for 
numerical calculation, where s = •\/47r\E', p = a/Itt^E'^u, q = V^'^,v, and A = . Their paper gives deflnitions and 
the full set of field equations for these quantities. The physical quantities used in our discussion can be expressed in 
terms of these quantities as follows: 





(Al) 






87ra2 ' 


(A2) 




(A3) 


pq 
47ra2' 


(A4) 


fq + gp 

a^la^pq 


(A5) 
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We use a two-step finite difference scfieme sliglitly different from tlie one proposed by Hamade and Stewart [28 
First we write tlie equations sclrematically as 



wlrere 



V,u = F{y,z), 
z,v = G(y,z), 



(A6) 
(A7) 



g 
f 



At the first step we predict the values at n{u, v) using those at v — h) and e(u — h, v) as 

Vn = ye + hF{ye,Ze), 



Zw + -jh{G{ytn,Zw) + G{yn, z„)). 



At the second step we can correct the prediction via 

1 



Vn = -{yn+ye + hF{yn,Zn)), 
Zn — -;^{Zn ~^ Z^ + hG{yn, Zn))- 



(A8) 
(A9) 

(AlO) 
(All) 



This is the scheme adopted by Due to Eq. (|A11|I . however, the scheme is only accurate to first-order, so we use 
instead 



Zn — Zn , 



(A12) 



since this is accurate to second-order and does not suffer from numerical instabilities. 

In general, Eq. (|A9|I can be implemented only implicitly. However, in this special case, it is possible to determine 
Zn explicitly as follows: 



"t" 2 ^Qn 5 



1 - hd, 



gn 




1 ^ /«<?« 

2 fn 



1 + 1/.^ 
2 fn 



(A13) 
(A14) 

(A15) 
(A16) 
(A17) 



where we put 



(A18) 



Note that the order of determination is very important here. 
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APPENDIX B: EXACT SOLUTIONS IN THE DOUBLE-NULL FORMULATION 



1. Flat Friedmann solution 



The flat Friedmann solution with a massless scalar field can be written as 

ds^ ^a^ [- dr]^ + dx^ +X^{dO^ + sin^ ddcj)^ )], (B 1 ) 
where — C^rj and C is a positive constant. This solution can be rewritten in double-null coordinates as 

ds^ = -a^dudv + r^ide^ +sirL^6d(t>^), (B2) 

a' -C^m^), (B3) 



I \ 1/2 
U + V ^ 



^L.f'^V.o, (B5) 



2 V 2 , 

where Sq is an arbitrary constant, and u and v are related to rj and x through 

u^VX, v = T] + x, (B6) 
respectively. The cosmological time t, where dt — adrj, t is given by 

t^\cr,'''^^^C{u + vf'\ (B7) 



1 / I \ -3/2 
1 U + 



The Hubble parameter H — d{hi a) /dt is given by 

There are two important horizons in this spacetime: the cosmological particle horizon and the cosmological apparent 
horizon. In the above coordinates, the cosmological particle horizon is given by u = or ry = while the cosmological 
apparent horizon is given by 3u + f = or 77 = x/2- Hence, the cosmological apparent horizon is spacelike and outside 
the cosmological particle horizon. The conformal diagram of the flat Friedmann solution is described in Fig. Q 

The physical quantities are 



C jv- u)^ 
16V2 (w + m)3/2 ' 



(B9) 



3 1 

= FT^T — \ — ^3' (^^0) 
J = 0, (Bll) 



87rC2 {u + v)3 ' 

1 V — u 

2V2V + U' 



1 V — u , , 

V = . (B13) 



This solution admits the following initial data: 



siuo,v) = ylln^^)+so, (B14) 
s{u,vo) = ^In (^ "^^° ) +so, (B15) 



2 V 2 



a'iuo,v) = (B16) 
aHu,vo) ^ C'(^), (B17) 
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and 



r{uo,vo) 
f{uo,vo) 
9{uo,vo) 



C 



uo + Wo 



r{uo,vo) 
r{uo,vo) 



1/2 



1 



Vo - Uq 



2{uo+vo) vo-uo 

1 1 

+ 



2{uo + vo) vo-uo 



(B18) 
(B19) 
(B20) 



2. Schwarzschild solution 



The Schwarzschild solution is given by 



2M\ ^ 



This can be rewritten in double-null coordinates as 

2M\ 



1 - -^j dudv + r^{d0^ + sin^ edcj)'^), 



where u and v are given by 



dn f, 2M\-' 

u = t-r^, V = t + r*, —— =1 

ar \ r J 



(B21) 



(B22) 



(B23) 



However, as is well known, this coordinate system has a coordinate singularity on the black hole event horizon. We 
therefore consider another coordinate system. For example, we can prescribe the initial data for this solution as 
follows: 

(B24) 
(B25) 
(B26) 

(B27) 
(B28) 





= a'^{u,vo) 


= const 


s{uo,v) 


= s{u,vo) 


= So, 


r{uo,vo) 


= ro, 




f{uo,vo) 


1 

~ ~2' 




9{uo,vo) 




2M\ 
ro J ■ 
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TABLE I: Model parameters and initial mass within black hole event horizon compared to cosmological apparent horizon 



Models 


Mm 


Am 


Mo 


Ml 


Mo 


Ml 


J^tBHEH/niCAH 


A 


-1/3 


0.02 


-0.5 


1.1 


1 


4 


0.972 


B 





0.02 


-0.5 


1.1 


1 


2 


0.223 


C 


-1/3 


0.5 


-0.5 


1.1 


1 


2.5 


0.399 


D 


-0.5 


0.5 


-0.5 


1.1 


1 


3.5 


0.727 



TABLE II: Initial and final radii of black hole horizons and their ratios 



Models 


BHAH v = vo 


BHAH M = Ml 


Ratio 


BHEH M = Mo 


BHEH M = Ml 


Ratio 


A 


0.373 


0.668 


1.79 


0.377 


0.694 


1.84 


B 


0.0858 


0.124 


1.44 


0.0985 


0.125 


1.27 


C 


0.154 


0.232 


1.51 


0.173 


0.241 


1.39 


D 


0.279 


0.470 


1.68 


0.300 


0.485 


1.62 



TABLE III: Initial and final masses within black hole horizons and their ratios 



Models 


BHAH M = Mo 


BHAH M = Ml 


Ratio 


BHEH M = Mo 


BHEH M = Ml 


Ratio 


A 


0.187 


0.334 


1.79 


0.187 


0.334 


1.79 


B 


0.0429 


0.0614 


1.43 


0.0429 


0.0614 


1.43 


C 


0.0769 


0.116 


1.51 


0.0769 


0.116 


1.51 


D 


0.140 


0.235 


1.68 


0.140 


0.235 


1.68 
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initial singularity 



FIG. 1: The conformal diagram of the flat Priedmann spacetime with a massless scalar field. The cosmological apparent 
horizon is spacelike and outside the cosmological particle horizon. 



BH singularity 




initial singularity 

FIG. 2: The conformal diagram of the possible causal structure of the calculated model. The region enclosed by a diamond is 
the calculated region [wo,mi] x [uo,t'i]. 




FIG. 3: Schematics figure of the initial data setting. The relation between Um and CAH is true only for Models A-C. 
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V V 

(c) (d) 

FIG. 4: The locations of the black hole event horizon, black hole apparent horizon, and cosmological apparent horizon in the 
{u, v) plane for Models A-D are plotted in (a) (d), respectively. For Model B, the location of the cosmological particle horizon 
of the unperturbed flat Friedmann solution is also plotted. 
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FIG. 5: The evolution of the area radius r of the black hole event horizon, black hole apparent horizon and cosmological 
apparent horizon for Models A- D is plotted in (a) -(d). For Model B, the radius of the cosmological particle horizon of the 
unperturbed flat Friedmann solution is also plotted. 
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FIG. 6: The evolution of the mass rn contained within the black hole event horizon, black hole apparent horizon and cosmological 
apparent horizon for Models A D is plotted in (a) (d). The curves for the black hole event horizon and the black hole apparent 
horizon are indistinguishable. For Model B, the mass contained within the cosmological particle horizon of the unperturbed 
flat Priedmann solution is also plotted. 
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FIG. 8: The area radius r along (a) ingoing null rays v — const and (b) outgoing null rays u — const for Model A. The arrows 
at u = Um = —1/3 and u — ubheh — 0.702 in (a) denote the matching outgoing null surface and black hole event horizon, 
respectively. 
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FIG. 9: Same as Fig.lHbut for 2m/r. 
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FIG. 10: The snapshots of scalar field profiles for Model A on a sequence of hypersurfaces with (a) constant v and (b) constant 
T. The profiles are indicated by broken curves for the numerical simulations and by solid light curves for the flat Friedmann 
solution. The heavy solid curves give the position of the event horizon. 
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FIG. 11: The snapshots of the energy density pu and the momentum density J of the scalar field measured by the observer 
moving normal to the T = (u + v) /2 = const spacelike hypersurface for Model A. Note that the event horizon expands with 
respect to r as time proceeds. 



24 



0.5 



1.5 



Event horizon 
v=1.00 
v=2.00 
v=3.00 
v=4.00 



2.5 



(a) 



1.5 
1 

0.5 


-0.5 
-1 

-1.5 
-2 



Event horizon 
v=1.00 
v=2.00 
v=3.00 
v=4.00 



0.5 



1.5 



2.5 



(b) 



FIG. 12: The snapshots of (a) the energy density ps and (b) the velocity V of the stiff fluid equivalent to the scalar field for 
Model A. It is noted that both ps and V are observer-independent. The event horizon always expands as time proceeds. 



